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ABSTRACT 


The  excess  volume  function  for  the  benzene  ♦  trichloro¬ 
ethylene  (TCE)  system  was  studied  at  temperatures  of  10,  25, 
30,  and  40  °C  in  this  investigation.  Prior  to  describing 
this  system,  benzene  «■  cyclohexane  at  25  °C,  the  test 
system,  was  analyzed  to  determine  the  ability  of  the 
investigator  to  obtain  accurate  results  with  the  apparatus 
and  procedures  described  in  this  work.  The  results  of  the 
test  system  were  in  excellent  agreement  with  those  of  other 
experimenters  and  with  the  parameters  described  by  Handa  and 
Benson  in  their  summary  of  this  system. 

The  excess  volume  function  of  each  system  studied  was 
described  by  a  smoothing  equation  that  was  derived  from  an 
average  of  72  data  points  per  system.  All  of  these  systems 
were  found  to  have  excess  volume  curves  that  were  uniformly 
positive  in  shape.  The  benzene  TCE  system  at  all  tempera¬ 
tures  had  maximum  values  at  about  58  percent  benzene.  This 
maximum  value  for  the  excess  volume  was  related  to  tempera¬ 
ture  by  the  equation  V  =  -.00140  *  T  ♦  .653697,  where  T  is 

in  K  and  V  _  is  in  cc/mol.  The  data  were  also  checked  for 
eT 

consistency,  by  the  data  point,  by  the  data  set,  and  by  the 
system  at  a  given  temperature,  with  favorable  results.  The 
relationship  between  temperatures  and  excess  volumes  was 
noted  to  invert  near  both  ends  of  the  excess  volume  curve. 
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CHAPTER  II 


BACKGROUND 

For  the  purposes  of  this  discussion,  a  binary  mixture 
will  be  defined  as  a  physical,  non-reactive  combination  of 
two  components  in  a  single  fluid  phase.  There  exist  thermo¬ 
dynamic  properties  of  each  particular  mixture  that  distin¬ 
guish  it  from  any  other  mixture  of  the  same  two  components. 
These  properties  can  be  divided  into  two  broad  categories: 
intensive  and  extensive.  An  extensive  property  is  one  that 
depends  on  the  amount  of  material  present.  These  properties 
include  total  mass  <M),  total  volume  <V>,  enthalpy  (H), 
entropy  <S),  and  Gibbs  free  energy  (G).  An  intensive 
property,  on  the  other  hand,  is  independent  of  the  amount  of 
mixture  present.  Temperature  <T>,  total  pressure  (P), 
specific  volume  (v),  and  density  <^o )  are  examples  of 
intensive  properties.  Due  to  their  intensive  nature  and 
ease  of  measurement,  most  of  the  other  properties  are 
described  as  functions  of  temperature  and  pressure. 

Additionally,  although  much  Information  is  generally 
known  about  single  component  system  properties,  relatively 
little  ie  available  on  the  properties  of  mixtures. 

Therefore,  it  is  helpful  to  discuss  the  properties  of 
mixtures  in  terms  of  the  pure  component  properties.  One  way 
of  doing  this  is  to  define  a  reference,  or  standard,  state 
as  the  pure  component  at  the  temperature  and  pressure  of  the 


system 


This  standard  state  is  referred  to  as  "Raoults  Law 


standard  state  (Prausnitz,  pp.  21-22)  and  will  be  used 
throughout  this  work. 

For  an  ideal  solution  to  exist,  the  total  volume  of  the 
solution  must  be  equal  to  the  sum  of  the  total  component 
volumes.  In  terms  of  molar  volumes: 

yID  *  ZX  \°  ,2> 

where:  V  is  the  total  molar  volume  of  an  ideal 

solution 

V^°  is  the  molar  volume  of  component  i  in 
solution 

x^  is  the  mole  fraction  of  component  i  in 
solution 

Ideal  behavior  is  seldom  observed  in  real  solutions. 
Therefore,  a  correction  factor  needs  to  be  included.  (Smith, 
p.  245) 

V  =  ^x±  V±0  +  (3) 

where:  V  is  the  total  molar  volume  of  a  real  solution 

V  is  the  molar  volume  change  on  mixing 
A  measure  of  the  deviation  of  a  solution  from  ideal 
behavior  is  the  excess  volume,  which  is  defined  as:  (Smith, 


266) 


where : 


.ID 


is  the  total  molar  excess  volume  of  a 
solution 


( 4  > 


It  is  convenient  at  this  time  to  introduce  the  concept 

of  partial  molar  volumes.  (Prausnitz,  p.  196) 

v\  =  h  nV/  9  n,  I  ,,  _  .  (5) 

i  I  il  T,P,nj 

:  V  is  the  partial  molar  volume  of  component  i 


where 


Note:  This  volume  is  calculated  at  conditions  of 


constant  temperature,  pressure,  and  total  moles  of  all 


components  except  that  of  component  i. 


This  concept  is  particularly  useful  in  that  for  all 


solutions  at  constant  temperature  and  pressure: 


v  -2-i  vi 


Now,  referring  to  equations  3  and  4, 


-  ‘ZVi0  -4v>  -£*^1° 


Further  simplification  yeilds 


VE  =  £  V 


Equation  8  results  since,  by  definition.  A  =  0, 


(Prausnitz,  p.  195)  and  thus,  the  excess  volume  is  not  a 


'new"  property.  Other  excess  properties,  for  example,  the 


excess  Gibbs  free  energy,  are  different  from  the  change  in 


the  property  on  mixing. 


By  definition  of  Raoult's  standard  state, 


V  ID  =  V  ° 
i  i 


Since  equation  S  applies  to  all  mixtures  and  with 


equations  2  and  9, 


yID  mE*±*±  "Xx±w±t 


is  another  way  of  defining  an  ideal  solution. 


E  A 

Since  V  * £1V,  if  the  molar  volume  change  on  mixing 


can  be  determined  experimentally,  so  can  the  excess  molar 


volume  of  mixing.  Thus,  a  measure  of  the  deviation  of  a 


solution  from  ideal  behavior  can  be  determined. 


Experimental  data  for  molar  excess  volumes  of  binary 


KM 


solutions  are  usually  fit  to  a  smoothing  function  of  form: 


(Handa,  p.  2021) 


*1  (1 


x€  >  Z! 


*  j»0  3 


a  ,  <  1  -  2x.)J  (11) 


where:  a  are  the  fitting  constants  determined  by  a 
J  least-squares  analysis 


K 

is 


n  is  the  number  of  coefficients 


The  units  of  V  are  generelly  cc/mol. 


The  number  of  coefficients  to  be  used  in  smoothing  the 
data  is  determined  by  carrying  out  the  least-squares 
analysis  for  increasingly  larger  values  of  n.  Usually,  n 
varies  from  2  to  5.  For  each  fit,  a  standard  deviation  is 
calculated. 


n 

rr  =  t  Z*  <  v . 


Z  (V.E  observed  -  V  E  smoothed)2  / ( m  -  n)]i/2  (12) 

H  i  i 


where:  m  is  the  number  of  experimental  data  points 

The  minimum  number  of  coefficients  is  selected  for  the 
smoothing  equation  that  will  yield  the  lowest  population 
standard  deviation. 

This  curve  fitting  technique  is  extremely  flexible  and 

capable  of  describing  most  of  the  possible  excess  molar 

volume  curve  types.  As  noted  in  Figures  1,  2,  and  3,  excess 

volume  data  can  be  uniformly  positive,  uniformly  negative, 

or  can  have  both  b  minimum  and  a  maximum  while  changing  sign 

somewhere  in  the  middle  of  the  mole  fraction  range. 

Although  quite  a  large  number  of  binary  systems  have 

been  studied,  relatively  few  oi  them  have  been  investigated 

over  a  wide  range  of  temperatures.  The  bulk  of  the  excess 

o 

volume  research  has  been  conducted  near  25  or  30  C.  For 


m 

m 

H 

BS 


w 


those  systems  that  were  studied  at  several  temperatures, 
there  appears  to  be  no  discernable  trend  of  the  data  with 
increasing  temperature. 

The  purpose  of  this  thesis  was  twofold.  First,  to 
investigate  the  benzene  +  cyclohexane  system  at  25  °C  for 
the  purpose  of  establishing  a  comparable  level  of  accuracy 
for  the  procedures  and  equipment  used  in  this  investi¬ 
gation.  And  second,  to  investigate  the  benzene  +  trichloro¬ 
ethylene  system  at  a  variety  of  temperatures  to  determine 
the  significance  of  temperature  change  on  this  previously 
unstudied  system. 

As  an  aside,  this  investigation  also  focused  on 
improving  the  apparatus  and  procedures  previously  set  out  in 
an  effort  to  improve  the  accuracy  of  the  experimental  data. 
Where  appropriate,  the  significant  improvements  and  lessons 
learned  will  be  noted. 
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CHAPTER  III 


INSTRUMENTATION  AND  EQUIPMENT 

The  dilatometer  used  in  this  invest igat ion,  shown  in 
Figure  4,  was  essentially  that  designed  by  Kumeran  and 
McGlashan  (Kumaran,  pp.  260-262),  with  a  few  modifications. 
The  capillary  tubes,  noted  as  C,  Cl,  and  C2,  were  precision 
bore  capillaries  with  a  . 05  millimeter  inside  diameter. 
Capillary  C  had  a  reference  mark  CR  near  the  bottom.  The 
buret,  labeled  B,  was  one  centimeter  precision  bore  tubing. 
The  bulb  at  the  bottom  of  the  buret,  BB,  was  a  blown  glass 
bulb  with  a  volume  of  about  10  cc  below  the  reference  mark 
BR.  Mixing  bulb  A  and  tubes  SI  and  S2  were  not  calibrated. 
The  taps,  T1  and  T2,  were  Teflon  and  had  Buna-N  O-rings 
midway  up  the  shaft.  Tap  T2  also  had  a  level  bubble 
attached  to  its  screw  cap.  The  small  addition  on  top  of 
capillary  C  was  half  of  a  ground  glass  joint  that  was 
attached  to  capillary  C  by  a  small  piece  of  Tygon  tubing. 

The  back-pressure  device  discussed  later  was  attached  to  the 
dilatometer  at  this  joint.  The  dilatometer  was  assembled  by 
Custom  Glassblowing  of  Louisville. 

A  plexiglass  frame  shown  in  figure  5  was  built  to 
support  the  dilatometer.  Two  of  the  three  support  blocks  on 
the  frame  were  adjustable  from  the  top  so  the  dilatometer 
could  be  leveled  while  it  hung  in  the  temperature  bath.  A 
thermometer  was  also  affixed  to  the  stand  for  monitoring  the 
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resistance- type  element  thBt  was  used  only  as  a  boost  heater 


to  help  bring  the  bath  initially  to  a  temperature  near  the 
set  point.  There  was  also  a  propeller  and  baffle  assembly 
to  circulate  the  water  and  keep  the  temperature  uniform 
throughout  the  bath.  The  temperature  was  controlled  by 
means  of  a  thermoregulator  which  protruded  through  the  top 
of  the  bath.  By  carefully  setting  the  controls  on  the  bBth, 
one  could  control  the  temperature  to  +.  .  001  °C  of  the  set 
point  value. 

When  cooling  was  required,  the  Neslab  PBC-2  II  portable 
bath  coolerin  Figure  8  was  used  in  conjunction  with  a 
Cryotrol  temperature  controller.  The  cooling  probe  was 
inserted  through  an  opening  in  the  top  of  the  temperature 
bath  and  the  temperature  sensor  was  inserted  through  a  small 
hole  in  the  top  of  the  dilatometer  support  stand. 

All  weights  were  measured  with  an  Ainsworth  model  Delta 
N-V  (Type  24N )  scale.  This  scale  had  a  digital  readout  to 
♦.  .  00001  grams,  with  the  last  two  digits  being  read  on  a 
vernier  scale. 

Height  measurements  were  made  with  the  use  of  cathetom- 
eters.  A  Gaertner  Scientific  Corporation  precision  cathetom 
eter  was  used  to  measure  all  heights  except  those  on  the 
buret.  This  cathetometer  could  be  read  to  ♦  .0001  cm,  with 
the  last  digit  being  read  on  a  vernier  scale.  The  height 
measurements  on  the  buret  were  made  with  an  Eberbach  cathe- 
tometer  that  could  only  be  read  to  ♦  .01  cm,  with  the  last 
digit  being  read  on  a  vernier  scale. 
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required  to  reach  thermal  equilibrium. 

Next,  with  the  buret  on  the  right,  the  dilatometer  was 
tilted  clockwise  such  that  a  small  amount  of  mercury  spilled 
into  the  buret.  The  dilatometer  was  then  placed  back  in 
the  temperature  bath  and  thermal  equilibrium  was  again 
reached.  All  mercury  levels  were  then  recorded  with  respect 
to  the  appropriate  reference  mark. 

A  pressure  correction  was  then  conducted  by  placing  a 
known  back-pressure  on  the  system  and  measuring  its  effect. 
This  correction  is  an  essential  factor  in  the  excess  volume 


data  for  two  reasons.  First,  there  may  be  some  small  amount 
of  gas  bubbles  entrained  in  the  system  which,  if  not 
corrected  for,  could  expand  and  influence  the  data.  Also, 


with  taps  T 1  and  T2  sealed,  pressure  differentials  can  be 
set  up  within  the  dilatometer  by  virtue  of  the  differing 
mercury  column  heights. 

The  tilting  of  the  dilatometer  and  measurement  of  the 
mercury  levels  was  repeated  successively  until  the  capacity 
of  the  dilatometer  was  reached.  Then  the  roles  of  the 
solvent  and  the  solute  were  reversed,  and  the  experiment 
was  repeated. 

Finally,  the  data  were  evaluated  for  internal  consis¬ 
tency  as  well  as  consistency  between  data  sets.  If  these 
evaluations  proved  unsatisfactory,  additional  data  were 
taken  as  required  to  achieve  the  desired  level  of  agreement. 


u Bed  for  the  benzene  were:  .8895  cc/mol  at  10  C,  .8790 
cc/iriul  at  25  °C,  .8884  cc/mol  at  30  °C,  and  .8577  cc/mol  at 
40  °C.  (API  44  Tables). 

The  trichloroethylene  used  here  was  a  99  ACS 

reagent  grade  product  of  Aldrich.  The  densities  used  here 
were:  1.4825  cc/rnol  at  10  °C,  1.4554  cc/mol  at  25  °C, 

1.4475  cc/mol  at  30  °C,  and  1.4300  cc/mol  at  40  °C. 

(Gallant,  p.  42). 

The  cyclohexane  used  in  this  work  was  a  99.  9  +  '/.  pure 
product  of  Burdick  and  Jackson.  At  25  °C,  the  density  of 
cyclohexane  is  . 7739  cc/mol. 


CHAPTER  V 


RESULTS 


A.  Data  Analysis. 

A  total  of  19  runs  w ere  made  for  this  thesis,  three  for 
the  test  system,  benzene  ♦  cyclohexane,  and  16  for  the 
benzene  ♦  trichloroethylene  system.  The  experimental  data 
for  each  run  were  initially  written  in  a  data  log.  The 
purpose  of  run  one  was  primarily  to  get  familiar  with  the 
operating  requirements  of  the  system.  The  data  for  this  run 
were  discarded  for  numerous  procedural  errors  that  obviously 
effected  the  results.  The  data  for  runs  number  9,  10,  and 

13  were  also  discarded  because  various  parts  of  the  dilatom- 
eter  broke  during  each  run.  In  run  9,  tap  T1  was  broken 
into  two  parts  under  the  O-ring.  In  run  10,  the  magnetic 
stirrer  initially  located  inside  the  mixing  bowl  developed  a 
crack  and  had  to  be  removed.  In  run  13,  capillary  C  broke 
where  it  joined  the  horizontal  tube.  Additionally,  runs  9 
and  12  were  discarded  because  of  degassing  problems  which 
left  large  gas  bubbles  inside  the  buret  after  the  run 
had  started.  No  further  data  analysis  was  performed  on 
these  five  runs. 

The  data  for  the  remaining  14  runs  were  transcribed  to 
a  computer  spread  sheet  for  further  analysis.  These  data 
are  located  in  Appendix  B.  With  the  aid  of  the  spread 
sheet,  the  mole  fraction  and  excess  volume  data  could  be 


calculated  immediately  after  each  data  point  was  recorded, 


as  shown  in  Appendix  C-i.  These  data  are  reported  here 
without  regard  for  significant  figures,  which  will  be 
addressed  separately.  Generally,  two  sets  of  data  were 
required  to  sufficiently  specify  the  excess  volume  curve  for 
a  given  temperature  and  set  of  components.  However,  late  in 
runs  15  and  18,  a  small  amount  of  solute  escaped  from  the 
mixing  bowl  during  a  tilt.  This  caused  a  separation  in  the 
mercury  column  in  capillary  C  and  would  also  have  resulted 
in  a  miscalculation  of  the  amount  of  solute  added  to  the 


mixing  bowl  if  the  runs  were  continued  beyond  that  point. 
This  error  apparently  resulted  from  rotating  the  dilatometer 
too  far  during  the  tilt.  Because  of  this  problem,  three 
sets  of  data  were  required  to  specify  the  curves  of  the 
benzene  «■  trichloroethylene  system  at  10  and  40  °C.  Run  7 
was  also  terminated  prematurely  due  to  degassing  problems, 
but  there  were  several  good  data  points  available,  so  runs  7 
and  8  were  combined  for  further  analysis. 

The  remaining  data  sets  were  grouped  by  temperature  and 


component  systems.  Set  one  was  benzene  ♦  cyclohexane  at  25 
o. 


C,  runs  2  and  3.  Set  two  was  benzene  +  trichloroethylene  at 


25  C,  initially  runs  4  and  5,  although  set  two  was  eventu¬ 


ally  comprised  of  runs  5  and  6.  Set  three  was  benzene 


trichloroethylene  at  30  C,  runs  7,  8,  and  11.  Set  four  was 


benzene  ♦  trichloroethylene  at  40  C,  runs  14,  15,  and  IS. 

o 


Finally,  set  five  was  benzene  ♦  trichloroethylene  at  10  C, 
runs  17,  18,  and  19.  The  mole  fraction  and  excess  volume 
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data  from  these  sets  were  converted  into  data  files  using 
the  program  "Makefile".  These  files  were  then  accessed  by 
the  program  "Curvfit"  which  yielded  the  curve  fitting 
constants  as  described  in  the  Background  section  of  this 
work.  A  listing  of  these  programs  is  found  in  Appendix  D. 

The  program  "Curvfit"  was  run  for  curve  fitting  constants 
ranging  from  2  to  5.  The  set  of  constants  yielding  the 
lowest  population  standard  deviation  was  then  used  to 
describe  the  smoothing  function.  These  data  are  presented 
in  Appendix  C-2. 

Additionally,  several  data  consistency  checks  were 
made.  First,  an  internal  consistency  check  was  made  using 
the  population  standard  deviation.  A  population  standard 
deviation  above  that  reported  by  other  investigators  meas¬ 
uring  excess  volume  data,  greater  than  .01,  would  warrant  a 
set  of  runs  being  discarded.  Next,  consistency  between  data 
runs  in  a  set  was  investigated.  If  the  data  points  in  the 
overlap  region  between  the  two  runs  were  offset,  the  more 
suspect  run  would  be  done  over,  and  a  new  pairing  made. 
Finally,  consistency  between  successive  data  points  was 
checked.  This  was  accomplished  by  eliminating  the  data 
point  that  differed  the  most  from  the  smoothing  function  and 
re-running  the  "Curvfit"  program.  Then,  the  next  worse  data 
point  would  be  removed,  and  so  on  until  ten  percent  of  the 
data  points  had  been  removed,  one  at  a  time.  If  a  clus¬ 
tering  of  bad  data  points  was  observed,  the  run  in  which 
this  clustering  occurred  would  be  thrown  out.  This  was  the 


The  worst  condition  number  reported  in  this  work  was  593  for 


the  benzene  trichloroethylene  system  at  40  °C.  This  would 
indicate  that  2  significant  digits  might  be  eliminated. 
However,  it  should  be  recalled  here  that  the  cathetometer 
used  to  read  the  height  measurements  on  the  buret  was  only 
accurate  to  within  .01  cm.  For  the  data  reported  here,  this 
level  of  accuracy  translates  to  only  three  significant 
digits  for  mole  fractions  of  benzene  lees  that  . 20  or 
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greater  than  . 80,  and  four  significant  digits  for  the  rest 
of  the  data  range.  This  indicates  that  the  ill  conditioning 
of  the  coefficient  matrix  did  not  affect  the  results  of  this 
investigation,  but  might  well  affect  the  results  of  further 
investigations  if  the  accuracy  of  this  cathetometer  is 
improved. 


It  should  be  noted  here  that  when  all  five  systems 
studied  in  this  work  are  compared,  the  population  standard 
deviation  generally  decreases  when  the  number  of  data  points 
taken  increases.  Also,  for  the  same  number  of  curve  fitting 


constants,  the  condition  number  is  smaller  for  a  larger 


number  of  data  points. 

Minor  temperature  fluctuations  appeared  to  have  little 
effect  on  the  results.  Although  the  thermometers  used  to 


monitor  each  system  varied  from  a  National  Bureau  of 


Standards  calibrated  thermometer  graduated  in  . 05  C  incre¬ 


ments  for  the  25  C  runs,  to  in  uncalibrated  thermometer 


graduated  in  . 1  °C  increments  for  the  10  and  40  °C  runs,  the 


standard  deviations  show  no  trend  with  degree  of  accuracy  of 
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the  thermometer. 


This  is  attributed  to  the  fact  that  the 


bath  was  able  to  control  the  temperature  consistently  better 
than  any  of  the  thermometers  could  measure  it.  The  actual 
temperature  at  which  the  run  was  conducted  may  be  question¬ 
able,  but  the  consistency  of  the  data  due  to  temperature 
fluctuation  was  excellent.  Runs  1  through  6  were  conducted 
at  better  than  25  °C  ♦  .01  °C,  (thickness  of  the  marks  on 
the  thermometer  was  the  determining  factor  here).  Runs  7 
through  11  were  conducted  at  30  °C  *  . 05  °C,  due  to  lack  of 
thermometer  calibration.  Runs  12  through  16  were  conducted 
at  40  °C  ♦  • 1  °C  and  runs  17  through  19  were  conducted  at  10 
°C  ♦  . 1  °C,  again  due  in  part  to  the  lack  of  NBS  calibration 
and  also  to  the  level  of  graduation  being  . 1  °C. 

B.  The  benzene  +  cyclohexane  system. 

The  benzene  *  cyclohexane  system  at  25  °C  has  been 
proposed  as  the  test  system  for  excess  volume  data.  Such  a 
great  volume  of  data  has  been  reported  on  this  system  that 
it  is  possible  to  measure  the  accuracy  of  the  equipment  and 
experimental  technique  used  in  an  investigation.  Once  an 
accuracy  level  has  been  established  for  this  system,  an 
investigator  may  ascribe  the  same  level  of  accuracy  to  all 
new  data  he  generates  with  the  same  apparatus  and  technique. 

Handa  and  Benson  went  to  great  lengths  to  document  the 
test  system  parameters  in  their  work.  (Handa,  pp.  200-205) 
Table  1  contains  a  comparison  between  the  data  generated  in 
this  investigation,  the  summary  of  Handa  and  Benson,  and  a 
summary  of  14  other  works.  There  were  106  data  points  taken 
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in  this  work  in  an  effort  to  accurately  describe  this 


system.  The  average  number  of  data  points  used  in  the  14 
other  works  was  about  34.  Handa  and  Benson  used  a  composite 
of  three  works  Judged  by  them  to  be  the  best  and  most  repre¬ 
sentative  data  available.  There  were  164  data  points  total 
among  these  three  data  sets.  Comparison  1  in  Table  1  was 
based  on  a  least -squares  curve  fit  with  all  164  data  points 
weighted  equally.  Comparison  2  in  Table  1  was  based  on  a 
least -squares  curve  fit  with  each  of  the  164  data  points 
being  assigned  a  weight  equal  to  the  inverse  of  the  squared 
standard  deviation  of  that  point  from  the  smoothed  curve  in 
the  work  from  which  it  was  taken. 

TABLE  1 


DATA  COMPARISON  FOR  THE  BENZENE  ♦  CYCLOHEXANE  SYSTEM 


Std  dev. 

Location 

WE  , 

V  at  max. 

VE  at  .5 

( cc/ mol ) 

of  max. 

( cc/mol > 

( cc/mol ) 

This  work 

. 00275 

.  485 

.  6529 

.  6523 

Handa  L  Benson 
Comparison  1 

. 00157 

.  491 

.  6511 

.  6509 

Hanaa  &  Benson 
Comparison  2 

. 00161 

.  491 

.  6520 

.  6518 

Range  of  Other 

. 00003  - 

.  485  - 

.6315  - 

.6312  - 

Works 

.  0050 

.  494 

.  6555 

.  6552 

Figure  9  shows  the  smoothed  excess  volume  curve 

r 

calculated ?in  this  work.  The  equation  for  this  curve  is: 
VE  =  x1  (1  -  x^  1 2. 6092  ♦  .15423(1  -  2x^ ) 

*  .12096(1  -  2xx)2  -  .14636(1  -  2x1>3] 

The  experimental  data  upon  which  this  curve  was  based  are 


superimposed  on  the  figure  in  an  effort  to  show  hov  well  the 
smoothed  curve  actually  describes  the  data.  There  is  no 
point  in  presenting  a  figure  comparing  this  work  to  the  data 
presented  by  Handa  and  Benson  because  the  curves  are  indis¬ 
tinguishable  on  any  reasonable  scale. 

C.  The  benzene  +  tricloroethy lene  system. 

It  should  be  noted  first,  that  although  the  accuracy  of 
these  data  cannot  be  shown  directly  by  comparison  with  known 
values,  it  can  be  inferred  from  the  accuracy  of  the  benzene 
♦  cyclohexane  system  discussed  above. 

Figures  10  through  13  are  the  smoothed  curves 
describing  the  benzene  trichloroethylene  system  at  10,  25, 
30,  and  40  °C,  respectively,  as  calculated  in  this  work  from 
the  equation: 

VE  =  x,  (1  -  x , ) [ a  +  a,  <1  -  2x, >  ♦  a_  (1  -  2x,)2 

1  1  O  1  12  1 

♦  a^  <1  -  2x^)2  ♦  a^  < 1  -  2x^>*3 
Table  2  contains  the  curve  fitting  constants  used  to  cal¬ 
culate  these  curves,  and  the  standard  deviation  from  the  fit 
for  each  temperature. 

All  four  of  these  curves  are  overlaid  in  Figure  14  for 
comparison.  There  was  little  utility  in  describing  the 
curves  at  both  25  and  30  °C,  as  there  was  very  little  change 
between  the  two  sets  of  data.  In  fact,  throughout  much  of 
the  range  of  these  two  curves,  they  are  sufficiently  close 
that  the  standard  deviations  of  each  data  set  could  cause 
these  sets  to  become  intertwined. 


FIGURE  1  1 

BENZENE  +  TRICHLOROETHYLENE  25  C 


FIGURE  13 


All  four  of  these  data  sets  have  reasonable  standard 

o 

deviations,  but  the  10  and  25  C  sets  had  notably  higher 

standard  deviations  than  those  for  the  30  and  40  °C  sets. 

See  Appendix  C-2  for  specific  information  about  standard 

deviations,  data  point  error,  and  condition  numbers  for  each 

temperature  in  this  system.  Finally,  the  data  appear  to  be 

sufficiently  good  for  information  at  other  temperatures 

between  10  and  40  °C  to  be  interpolated  with  reasonable 

accuracy.  An  empirical  equation  relating  the  height  of  the 

excess  volume  peak  to  the  temperature  of  the  system  is: 

V  _  =  -.00140  *  T  +  .653697 
el 

Where  the  units  of  V  _  are  cc/mol  and  T  is  in  K.  This 

eT 

equation  was  derived  from  a  linear  regression  of  the  peBk 
values  for  the  four  temperatures. 


CHAPTER  VI 


CONCLUSIONS 

The  data  for  the  benzene  *  cyclohexane  system  compare 
well  with  other  experimenters  in  terms  of  standard 
deviation,  maximum  value,  and  general  shape  of  the 
curve,  but  is  slightly  low  in  the  location  of  the 
maximum  point  on  the  x  axis,  although  it  is  within  the 
range  of  acceptibility . 

The  data  presented  for  the  benzene  ♦  trichloroethylene 
system  are  as  good  as  the  data  for  the  benzene  * 
cyclohexane  system,  have  standard  deviations  that  are 
within  acceptible  limits,  have  peak  values  near  58 
percent  benzene,  and  can  be  used  to  interpolate  data  at 
other  temperatures. 

For  the  benzene  trichloroethylene  system,  as  the  temp 
erature  increases,  the  number  of  fitting  constants 
required  to  describe  the  curve  generally  increases. 

For  the  majority  of  the  benzene  trichloroethylene 
system,  the  excess  volume  of  mixing  decreases  with 
increasing  temperature,  although  an  inversion  of  this 
relationship  takes  place  near  both  ends  of  the  excess 
volume  versus  mole  fraction  curve. 

As  more  data  points  are  taken  to  describe  a  given 
system,  the  standard  deviation  end  the  condition 


decrease. 


The  inaccuracy  of  reading  the  temperature  from  an 
uncalibrated  thermometer  had  no  effect  on  the  internal 


consistency  of  the  data.  This  implies  that  the  temper¬ 
ature  bath  was  controlled  to  a  degree  of  precision 
better  than  any  of  the  thermometers  were  capable  of 
reading . 

With  an  increase  in  the  precision  of  the  cathetometer 
that  measures  the  buret  data,  an  increase  of  precision 
in  the  excess  volume  data  and  the  standard  deviations 
would  follow.  The  excess  volume  precision  would  likely 
increase  decimal  point  for  decimal  point  with  the 
precision  increase  of  the  cathetometer,  with  modifi¬ 
cations  made  for  ill  conditioning  if  too  few  data 
points  are  collected. 


The  curve  fitting  constants  and  standard  deviations  for 
each  of  the  benzene  trichloroethylene  systems  are  as 


noted  in  Table  2. 


CHAPTER  VII 


RECOMMENDATIONS 


Obtain  a  second  high  precision  cathetometer  to  measure 
the  height  data  on  the  buret. 

Obtain  NBS  calibrated  thermometers  to  cover  a  wide  range 
of  temperatures,  perhaps  at  ten  degree  increments  from  0 
to  100  °C,  to  accurately  set  the  temperature  bath  at  the 
desired  set  point. 

Take  as  many  data  points  as  possible  in  further  investi¬ 
gations  to  lower  the  population  standard  deviation  and 
to  keep  the  ill  conditioning  of  the  coefficient  matrix 
from  negating  the  positive  effects  gained  from  the  use 
of  higher  precision  equipment. 

Continue  the  study  of  the  benzene  +  trichloroethylene 


system  at  extreme  temperatures,  such  as  -10,  0  and  60 

°C,  to  see  if  the  observed  trends  continue  to  manifest 
themselves. 

Obtain  and  calibrate  a  spare  set  of  capillaries  and 


burets  for  use  in  case  of  accidental  breakage. 
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Buret  Calibration 


Clean  the  dilatometer. 

Fill  a  syringe  to  a  reasonable  level  with  clean  mercury, 
less  than  half  full  is  acceptable.  The  rationale  for 
this  amount  is  that  if  the  weighing  will  be  done  in  an 
enclosed  scale,  the  entire  syringe  ( needle,  body,  and 
plunger)  must  fit  within  the  confines  of  the  scale  when 
the  scale  doors  are  shut.  Never  attempt  to  over  fill 
the  syringe  and  attempt  to  get  an  accurate  weight  with 
the  scale  doors  open  because  the  room  air  currents  are 
sufficient  to  cause  errors  in  the  weight. 

3.  Weigh  the  syringe  and  record  this  weight. 

4.  Tilt  the  dllatometer  such  that  the  buret  is  inclined 
slightly  above  horizontal.  Slowly  Inject  the  mercury 
into  the  buret  while  keeping  the  opening  of  the  needle 
flat  against  the  side  of  the  buret.  The  mercury  drops 
should  slowly  run  down  the  wall  of  the  buret  and  pool  up 
at  the  bottom. 

5.  When  all  mercury  possible  has  been  injected  into  the 
buret,  withdraw  the  syringe  and  reweigh  it.  Care  must 
be  taken  to  ensure  that  no  mercury  remaining  in  the 
syringe  is  allowed  to  escape.  Record  this  weight. 

6.  Repeat  steps  2-5  as  many  times  as  needed  to  ensure  the 
mercury  level  in  the  buret  is  above  the  reference  mark 
at  the  lower  end  of  the  buret.  Better  results  can  be 
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achieved  if  the  level  is  as  close  to  the  reference  line 


as  possible,  as  long  as  it  is  above  it. 

7.  Visually  inspect  the  mercury  in  the  buret  to  ensure  that 
no  air  bubbles  are  visible.  If  a  bubble  is  detected, 
tilt  the  dilatometer  in  a  manner  such  that  the  trapped 
air  bubble(s)  may  rejoin  the  air  above  the  mercury. 

8.  Replace  the  tap  on  the  buret  and  tighten  it.  For 
calibration  of  the  buret  only,  the  tap  with  the  level 
bubble  should  be  used  on  the  buret  and  the  plain  tap  on 
the  mixing  bowl. 

9.  Place  the  dilatometer  in  the  stand  and  place  the  stand 
in  the  temperature  bath.  Insert  the  cooling  unit  temper 
ature  sensor  into  the  bath  through  a  hole  in  the  top  of 
the  stand. 

10.  Level  the  dilatometer  by  using  the  level  bubble  and 
making  appropriate  corrections  with  the  height 
adjustment  rods. 

11.  Allow  the  system  to  reach  thermal  equilibrium  at  the 
desired  operating  temperature,  which  should  be  held 
constant  throughout  the  calibration. 

12.  Read  the  height  to  the  reference  line  and  the  height  to 
the  meniscus  using  a  cathetometer .  Record  these 


readings. 

13.  Fill  a  syringe  with  a  small  alloquit  of  mercury.  The 
smaller  the  amount,  the  better  the  calibration  will 
be.  One-half  cc  works  well. 

14.  Remove  the  temperature  sensor  from  the  bath  and 


Dilatometer  Cleaning  Procedures 


1.  The  dilatometer  should  be  cleaned  prior  to  each  loading. 

2.  Pour  nitric  acid  into  the  buret  in  sufficient  quantities 
to  fill  the  buret  and  the  mixing  bowl.  A  funnel  with  a 
short  stem  should  be  used. 

3.  Tilt  the  dilatometer  back  and  forth  in  such  a  manner 
that  the  nitric  acid  comes  in  contact  with  all  of  the 
inner  surface  of  the  dilatometer,  to  include  the 
capillaries. 

4.  Pour  out  the  nitric  acid  and  fill  as  much  of  the 
dilatometer  as  possible  with  distilled  water. 

5.  Rinse  the  dilatometer  as  in  step  3,  then  pour  out  water. 

6.  Add  ethyl  alcohol  to  the  dilatometer  as  in  step  2,  then 
rinse  as  in  step  3  and  pour  out  the  alcohol. 

7.  Execute  step  &  again  with  acetone  instead  of  alcohol. 

8.  When  the  acetone  has  been  poured  out ,  connect  a  vacuum 
pump  to  capillary  C.  Allow  the  vacuum  to  operate  until 
tubes  SI  and  S2  as  well  as  taps  T1  and  T2  are  dry. 

9.  Turn  off  the  vacuum  and  insert  both  taps  tightly.  Then 
turn  the  vacuum  back  on  and  let  it  operate  until  the 
dilatometer  is  completely  dry  on  the  inside. 

Note:  It  is  helpful  to  conduct  the  vacuum  operation  with 

the  dilatometer  in  an  inverted  position.  This  prevents 
acetone  from  collecting  in  the  bulb  at  the  bottom  of  the 
buret  and  speeds  the  drying  process  considerably. 
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Dilatometer  Loading  Procedures 


Refer  to  Figure  4  for  the  nomenclature  used  in  this  section. 

1.  First,  clean  the  dilatometer. 

2.  Inject  mercury  into  the  mixing  bowl  through  tap  T2  with 
a  long  needle  and  a  small  tube  attached  to  a  syringe. 
Care  should  be  taken  to  keep  the  tube  opening  just  above 
the  level  of  the  mercury  to  ensure  that  no  air  bubbles 
become  entrapped  in  the  mercury  pool. 

3.  As  the  horizontal  tube  fills  with  mercury,  inspect  peri¬ 
odically  for  the  formation  of  air  bubbles.  If  any 
bubbles  form,  tilt  the  dilatometer  appropriately  to 
allow  the  air  bubble(s)  to  escape. 

4.  Continue  to  fill  the  mixing  bowl  with  mercury  until 
about  two  thirds  of  the  upper  hemispherical  portion  of 
the  bowl  are  full. 

5.  Put  both  taps  in  place,  but  do  not  fully  seat  them. 

Allow  T2  to  be  open  enough  so  that  mercury  can  be 
Injected  through  S2  into  the  mixing  bowl. 

6.  Place  the  dilatometer  in  its  stand  and  lower  it  into  the 
temperature  bath  which  is  set  at  the  run  operating 
temperature. 

7.  When  thermal  equilibrium  has  been  reached,  remove  the 
dilatometer  from  the  temperature  bath  and  inject  just 
enough  mercury  into  the  mixing  bowl  through  S2  so  that, 
when  tap  T2  is  closed,  no  air  will  be  trapped  inside, 


ft 


into  £2  or  a  pressurizing  effect  when  T2  is  closed. 

Close  T2  and  place  the  dilatometer  back  in  the  temper¬ 
ature  bath. 

Agai  n,  allow  for  thermal  equilibrium  to  be  reached,  then 
measure  the  heights  of  the  mercury  columns  in  C,  Cl,  and 
C2  with  respect  to  CR. 

When  it  is  desired  to  operate  the  dilatometer  near  or 
above  the  ambient  temperature,  both  the  solvent  and  the 
solute  should  be  degassed  prior  to  loading  them  into 
the  dilatometer.  This  can  be  most  conveniently 
accomplished  by  bringing  each  component  to  boiling  just 
prior  to  injecting  it  into  the  dilatometer. 

Remove  the  dilatometer  from  the  bath  and  open  both 
taps,  being  careful  not  to  allow  any  water  to  enter  the 
dilatometer.  Inject  the  solute  component  into  the 
buret.  This  is  accomplished  with  the  help  of  a  long 
needle  and  a  syringe.  Fill  the  buret  with  this  com¬ 
ponent  by  placing  the  opening  of  the  needle  against  the 
side  of  the  buret  and  slowly  injecting  the  solute  in 
such  a  manner  that  it  flows  down  the  walls  of  the  buret 
in  a  film. 

Continue  to  fill  the  buret  until  the  level  of  the 
solute  is  above  the  upper  openings  of  both  Cl  and  C2. 
Momentarily,  tilt  the  dilatometer  toward  the  buret  such 
that  mercury  drops  flow  from  both  Cl  and  C2.  When  the 
air  trapped  inside  C2  has  escaped  into  the  buret,  tilt 


the  buret  upright  again. 
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Tilt  the  dilatometer  in  such  a  manner  that  mercury 


flows  into  the  buret  only  from  Cl.  Allow  mercury  to 


flow  from  Cl  into  the  buret  until  the  bulb  at  the 


bottom  is  full  of  mercury  and  the  meniscus  of  the 


mercury  is  slightly  above  reference  line  BR. 


As  mercury  falls  into  the  buret,  the  solute  level  will 


rise.  The  excess  solute  should  be  removed  with  a 


syringe  through  tap  T1  from  time  to  time. 


Once  this  portion  of  the  procedure  is  complete,  ensure 


the  solute  level  is  even  with  the  lower  opening  of  SI. 


Replace  tap  T 1,  but  leave  it  slightly  vented  to  SI. 


Fill  the  vacant  space  in  the  mixing  bowl  with  the 


solvent  as  in  step  2.  Bring  the  level  of  the  solvent 


even  with  the  lower  opening  of  S2. 


Replace  tap  T2,  but  leave  it  slightly  vented  to  S2. 


Place  the  dilatometer  back  into  its  stand  and  lower  it 


into  the  temperature  bath. 


When  thermal  equilibrium  has  been  reached,  seal  both 


taps  and  measure  the  heights  of  mercury  in  columns  C, 


Cl,  C2,  and  the  mixing  bowl  with  respect  to  CR  and  the 


height  of  mercury  in  the  buret  with  respect  to  BR. 


If  any  of  the  capillary  readings  cannot  be  made  conve¬ 


niently,  their  levels  may  be  adjusted  in  the  following 


manner.  Remove  the  dilatometer  from  the  temperature 


bath.  Open  tap  T1  partially,  and  either  allow  more 


solute  into  the  buret  or  remove  a  few  drops,  as 


m 
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c.  Set  the  heating  capacity  knob  to  ate  lowest  value. 

d.  Allow  the  bath  to  operate  and  observe  the  temper¬ 
ature. 

e.  Adjust  the  heating  capacity  upwards  until  the  temp¬ 
erature  an  the  hBth  begins  to  climb. 

f .  Verify  that  the  bath  can  reach  the  desired  operating 
temperature  by  adjusting  the  set  point  upwards  so 
that  the  bath  temperature  climbs  above  the  desired 
operating  temperature.  If  the  temperature  quite 
climbing  before  the  operating  temperature  is 
exceeded,  continue  adjusting  the  heating  capacity 
upwards  until  this  goal  is  achieved. 

Set  the  temperature  controller  on  the  bath  lower  than 
the  desired  operating  temperature  and  allow  the  bath  to 
stabilize  there. 

Adjust  the  temperature  controller  very  slightly  upwards 
and  again  allow  the  bath  to  stabilize  in  temperature. 
Repeat  the  procedure  in  step  6  until  the  operating  temp¬ 
erature  is  reached. 

The  temperature  controls  are  now  set  for  a  run  and 
should  not  need  adjusting  unless  ambient  conditions 
change  significantly. 

At  the  end  of  each  day  of  operation,  refill  the  temp¬ 
erature  bath  with  distilled  water. 

Let  the  bath  run  continuously  during  a  run  to  prevent 
instability  and  to  facilitate  the  next  day's  operation. 


Operating  Procedures 
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Before  beginning,  choose  an  operating  temperature  and 

set-up  the  temperature  bath  appropriately.  Also,  choose  the 

two  component  system  to  be  investigated. 

1.  Designate  one  component  as  the  solute  and  the  other  as 
the  solvent. 

2.  Clean  and  load  the  dilatometer,  recording  the  appro¬ 
priate  data  as  described  in  the  loading  procedures. 

3.  Remove  the  dilatometer  from  the  bath. 

4.  With  the  buret  on  the  right,  tilt  the  dilatometer  in  a 
clockwise  direction,  allowing  a  few  drops  of  mercury  to 
spill  into  the  buret  from  capillary  Cl.  Right  the  dila- 
tometer  and  place  it  back  in  the  temperature  bath. 

Note:  After  the  tilt,  inspect  capillary  C2  to  ensure  that 

there  is  no  separation  in  the  mercury  column.  If  this 
should  occur,  immediately  retilt  the  dilatometer.  An 
abrupt  tilt  reduces  the  risk  of  mercury  separation. 

5.  Allow  the  system  to  reach  thermal  equilibrium.  This 
takes  about  30  minutes,  and  can  be  visually  determined 
by  observing  the  mercury  in  capillary  C.  When  this 
mercury  level  reaches  a  stationary  location,  thermal 
equilibrium  has  been  reached. 

£..  Connect  capillary  C  to  the  back  pressure  device  and  turn 
the  stopcock  to  the  vent  position. 

7.  Read  and  record  the  levels  of  mercury  in  capillaries  C, 


and  C2;  the  mixing  bowl;  and  the  buret. 


Also, 


record  the  level  of  reference  marks  CR  and  BR. 

Turn  the  stopcock  to  the  closed  position  and  put  a  back 
pressure  of  about  24  cm  on  the  system  as  determined  by 
observing  the  change  in  height  of  the  mercury  in  the  U 
tube. 

Record  the  net  change  of  height  of  mercury  in  the  U 
tube,  ensuring  that  any  initial  differential  in  mercury 


column  heights  is  accounted  for. 


Record  the  change  in  height  of  mercury  in  capillary  C. 
For  the  first  few  tilts,  only  a  few  drops  of  mercury 
should  be  allowed  to  spill  into  the  buret.  This  facili¬ 
tates  a  more  accurate  determination  of  infinite  dilution 
values. 

After  the  first  few  tilts,  gradually  increasing  amounts 
of  mercury  may  be  spilled  into  the  buret  until,  even¬ 
tually,  the  change  in  mercury  height  in  the  buret  is 
about  1  cm  per  tilt.  Once  this  value  is  reached,  the 
run  should  be  finished  out  at  about  1  cm  of  mercury 
change  in  the  buret  per  tilt. 

This  procedure  should  continue  until  the  mercury  in  the 
buret  reaches  the  point  where  the  buret  cross-sectional 
area  begins  to  vary  as  it  merges  with  the  horizontal 
tube. 

Repeat  steps  1-13  with  the  solute  and  solvent  roles 
reversed.  This  procedure  should  yeild  two  curves  that 
merge  at  the  center  with  a  significant  overlap. 
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NOMENCLATURE  FOR  APPENDIX  B 


Mixing  bovl  on  the  dilatometer 
Buret  on  the  dilatometer 
Reference  mark  on  the  buret 
C 2  Capillary  tubes  on  the  dilatometer 

Reference  mark  on  capillary  C 
pi  Pressure  corrected  change  in  height  of  mercury 

in  C 

Back-pressure  applied  to  C 


EXPERIMENTAL  DATA  FOR  DILUTION  RUN  #2 
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CALCULATED  OflTfl  FOR  DILUTION  RUN  #19 


How  many  curve  fitting  constants  do  you  want? 
(must  be  between  2  and  5) 

4 


What  Is  the  name  of  the  data  file  to  be  used? 
bzcyc  1  o 

(Data  must  be  stored  in  x.y  pairs  with  format  f 7. 6 . lx , f 7 .6 . 
The  first  line  of  the  file  must  contain  the  number 
of  data  points  in  the  file,  with  format  13) 


gamma  1  = 
gamma  2  = 
gamma  3  = 
gamma  4  = 

x  benzene 


0 . 982598D+00 
0 . 9737250+00 
0.95819601-00 
0 . 9524790-100 
0 . 943940D+00 
0.8642180+00 
0.8567020+00 
0 . 844607D+00 
0.8209360+00 
0.8153880+00 
0.8034990+00 
0.7477670+00 
0 . 7365570+00 
0.7237920+00 
0.714  135D+00 
0 . 708  145D+00 
0.6811880+00 
0.6740670+00 
0.6543240+00 
0.6480910+00 
0 . 6447 10D+00 
0.6398170+00 
0 . 636642D  +  00 
0 . 6246070+00 
0.6188490+00 
0.6112730+00 
0 . 6051350+00 
0.5982170+00 
0.5908360+00 
0.5851490+00 
0.5795590+00 
0.5704540+00 
0.5645600+00 
0.5564890+00 
0.5510400+00 
0.5437760+00 


2 . 60920 
0 . 15423 
0 . 12096 
-0 . 14836 


mo  I  a  r  Ve 


0.4898500-01 
0.6833200-01 
0 . 1088800+00 
0. 1168510  +  00 
0. 1394680  +  00 
0.3038520+00 
0.3193630+00 
0.3405570+00 
0.3853290+00 
0.3853320+00 
0.4095990+00 
0.4916650+00 
0.5035080+00 
0.5172390+00 
0.5316830+00 
0.5320710+00 
0.5642260+00 
0.5656010*00 
0.5837540+00 
0.5852930+00 
0.5887410+00 
0.5956410+00 
0 . 5984740  +  00 
0.6066600+00 
0.6076570+00 
0.6138350+00 
0 . 613927D  +  00 
0 . 6224240+00 
0 . 6258330+00 
0.6270720+00 
0.6301770+00 
0 . 6380500+00 
0.6329180+00 
0.6422130+00 
0.6418090+00 
0.6462970+00 


0.462777D-01 
0 . 690231D-01 
0 . 1Q7496D+00 
0 . 121 243D+00 
0. 14 1367D+00 
0.3072550+00 
0.3209780+00 
0.3424 12D+00 
0.3820940+00 
0 . 390970D+00 
0 . 409457D+00 
0.4867 18D+00 
0 . 500435D+00 
0.5153280+00 
0 . 526083D  +  00 
0 . 532534D+00 
0.5594880+00 
0.5660430+00 
0 . 582986D+00 
0 . 587959D+00 
0 . 59058 1D+00 
0 . 594282D+00 
0 . 5966240+00 
0 . 605075D+00 
0 . 6088800+00 
0 . 613650D+00 
0 . 617318D+0Q 
0.6212410+00 
0 . 625 178D+00 
0 . 628036D+00 
0 . 630697D+00 
0 . 6347 14D+00 
0 . 637  104D+00 
0 . 640 108D+00 
0 . 64 1959D+00 
0 . 644 205D+00 


0 . 2707320-02 
-.6911160-03 
0. 1383870-02 

-  .4391700-02 
-. 1898590-02 
-.3402860-02 
-.1615280-02 
-.  1855170-02 
0.3235250-02 

-  .5637510-02 
0. 1420030-03 
0.4947170-02 
0.3072B1D-02 
0. 1911330-02 
0 . 560028D-02 
-.4625110-03 
0.473795D-02 
“ . 442270D-03 
0.767930D-03 
- . 266625D-02 
-. 1840420-02 
0. 135B79D-02 
0. 18499 lD-02 
0. 158466D-02 

-  .  1223120-02 
0. 1B4562D-03 

-  .339148D-02 
0. 1183020-02 
0.654957D-03 

-  .  96445  ID-03 
- . 520206D-03 
0.3335990-02 
-.4186040-02 
0 .2105230-02 

-  .  1501800-03 
0. 209188D-02 


Wm 


p: 


0.537856D+00 
0 . 52 18  6  ID*  00 
0.5152  140+00 
0.5053340+00 
0 . 4962490+00 
0 . 4898360+00 
0.4827310+00 
0 . 4759970+00 
0 . 172  1900-0  1 
0 . 3615800-0  1 
0 . 4946600-0  1 
0 . 7478100-01 
0 .8773600-01 
0 . 104eC3D+0Q 
0 . 1203580  +  00 
0  .  1303820+00 
0  .  1797180+00 
0 . 1909220  +  00 
0.2117270+00 
0.2247160+00 
0 . 230  720+00 
0 . 2422**. 0  +  00 
0 . 2509650+00 
0 . 2604000+00 
0.2721420+00 
0 . 280297D  +  00 
0 . 289399D+00 
0 . 2990450+00 
0.3082150+00 
0.3165610+00 
0 .3248880+00 
0.3309290+00 
0 . 33949  10  +  00 
0.3443720+00 
0 . 3528650+00 
0 . 3833  100  +  00 
0 . 3888590+00 
0 . 3968210+00 
0 . 4045940+00 
0.4124350+00 
0. 4187310+00 
0 .4262270+00 
0.4331450+00 
0 . 4397040+00 
0.4492720+00 
0 . 4538950+00 
0.4610030+00 
0.4665600+00 
0 . 473266D+00 
0.4776410+00 
0.4831680+00 
0.4904650+00 


0 . 6488830+00 
0 . 6497570+00 
0 . 6526760  +  00 
0 . 648356D  +  00 
0 . 65  10220  +  00 
0 . 6489330+00 
0 . 6533090  +  00 
0 . 65043 10+00 
0.4719900-01 
0.9593500-01 
0. 1258530+00 
0. 1836390  +  00 
0 . 2203940  +  00 
0 . 2620440+00 
0 . 2873660+00 
0.3130170+00 
0.4000600+00 
0.4189120+00 
0.4527780+00 
0 .4747  130  +  00 
0.4828850+00 
0.4934530+00 
0.508297D+00 
0 .5203160  +  00 
0.5377750+00 
0 . 5386730+00 
0 .549923D+00 
0 . 5604140+00 
0.5684270+00 
0.5735860+00 
0.5846160+00 
0.5879670+00 
0.5947550+00 
0 . 6006670+00 
0 . 6Q4277D+00 
0 . 6224450+00 
0.6249590+00 
0.6368910+00 
0 . 6377740+00 
0.6419540+00 
0.6397170+00 
0 . 645  108D  +  00 
0 . 6477730+00 
0 . 6480370+00 
0 . 6509320  +  00 
0 . 648650D+00 
0.6523110+00 
0 . 6535280  +  00 
0 . 6531480+00 
0.6516830+00 
0 . 6569130+00 
0 . 653824D+00 


0 . 645847D+00 
0 . 6494320  +  00 
0 . 650553D+00 
0.6518 18D+00 
0 . 652555D+00 
0 . 6528260+00 
0 . 652887D+00 
0 . 6527090+00 
0 .4632250-01 
0 . 954  1830-0 
0. 1287310+00 
0 . 18934  1D+00 
0 .2189400+00 
0 . 256447D  +  00 
0. 289146D+00 
0.3094640+00 
0 . 4007790+00 
0 . 4  194980  +  00 
0 . 4522780  +  00 
0 .4714390+00 
0 . 4800310  +  00 
0 . 495720D+00 
0 . 507 1 160+00 
0.5189500+00 
0.5329490+00 
0.542  199D  +  00 
0. 5520670+00 
0.5620030+00 
0 . 570953D+00 
0 . 578680D+00 
0 . 585995D+00 
0 . 59 1055D+00 
0.5978750+00 
0 . 60  1579D+00 
0 . 607707D+00 
0 . 626392D+00 
0 . 62925 1D+00 
0 . 63306 1D+00 
0 . 63645  2D  +  00 
0.6395450+00 
0 . 64  1792D+00 
0 . 644  192D+00 
0 . 646145D+00 
0 . 6477 66D+  00 
0 . 649 7 27D+00 
0 . 650505D+00 
0 . 65 1487D+00 
0 . 652073D+00 
0.6525720+00 
0 . 6527740+00 
0 . 652890D+00 
0 . 652809D+00 


0 . 3035950-02 
0. 325362D-03 
0 . 212270D-02 
-.3462320-02 
-. 1532550-02 

-  .  3893430-02 
0.4221910-03 

-  .2278180-02 
0 .8764860-03 
0  5166770-02 
- . 2877700-02 

-  .  5702400-02 
0 . 1454270-02 
0 . 559743D-02 
- .  1780200-02 
0 . 3553290-02 

-  .  71B890D-03 
-.5863470-03 
0.5003040-03 
0 .3273920-02 
0 . 2853730-02 
- . 226725D-02 
0. 118064D-02 
0 . 136638D-02 
0.4826190-02 
- . 352568D-02 
- .2144050-02 

-  .  15B892D-02 
-.  252556D-02 
-.509393D-02 
- . 1378620-02 
- . 308834D-02 
-.312049D-02 

-  .9122720-03 
-.3429540-02 
- . 394675D-02 
-.4291660-02 
0 . 382952D-02 
0. 132161D-02 
0.240685D-02 
- . 207453D-02 
0.9756750-03 
0. 1627610-02 
0 . 27  1440D-03 
0. 120475D-02 
- .  185520D-02 
0.824420D-03 
0. 1454620-02 
0.5762770-03 
- . 1090700-02 
0.4022570-02 
0 . 1015140-02 


*  benzene 

0.4957240+00 
0 . 6048330+00 
0.5102690+00 
0 . 6  164020+00 
0 . 5221560+00 
0.5291170+00 
0 . 5345800+00 


molar  Ve 

0 . 656270D+00  ' 
0.651586D+00 
0 . 654  1220+00 
0 . 6504950+00 
0.6538960+00 
0 . 6484580+00 
0 . 6468470+00 


0 . 6525850+00 
0 . 651870D+00 
0 . 651247D+00 
0 . 6503690+00 
0 . 6493770+00 
0 . 64  7 9  60D  +  00 
0 . 6466830  +  00 


0 . 368547D-02 
-  .  2836020-03 
0 . 2875270-02 
0 . 1262350-03 
0.451917D-02 
0 . 497927D-03 
0 . 164340D-03 


The  standard  deviation  for  this  set  of  fitting  constants  is:  0.2750330-02 
The  condition  number  for  the  A  matrix  is:  85.648 
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kv 
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How  many  curve  fitting  constants  do  you  want? 
(must  be  between  2  and  5) 


What  is  the  name  of  the  data  file  to  be  used? 


(Data  must  be  stored  in  x.y  pairs  with  format  f 7 . 6 , lx , f 7 . 6 . 
The  first  line  of  the  file  must  contain  the  number 
of  data  points  in  the  file,  with  format  13) 


gamma  1  = 
gamma  2  = 
gamma  3  = 


1.01164 

-0.41103 

0.07189 


x  benzene 


molar  Ve 


0 . 9690980*00 
0.9441360+00 
0.896863D+00 
0.84217  60+00 
0 . 785282D+00 
0 . 75 14 1 2D+00 
0. 705836D+00 
0.6691710+00 
0.6403360+00 
0.6077830+00 
0.6847370+00 
0.5590760+00 
0.5339840+00 
0.5247150+00 
0 . 5 15808D+00 
0.497848D+00 
0 . 487898D+00 
0 . 470640D+00 
0.4670300+00 
0.4524330+00 
0.4458520+00 
0 . 426399D+00 
0.4108870+00 
0.3763330+00 
0.3728130+00 
0 . 334998D+00 
0.3006540+00 
0.2655230+00 
0.2356530+00 
0.2002900+00 


0.4633100-01 
0.7992800-01 
0. 1247850+00 
0. 1679420+00 
0. 196586D+00 
0.2255080+00 
0 . 238624D+00 
0. 251752D+00 
0. 262429D+00 
0.2808830+00 
0.2767510+00 
0.2667200+00 
0.2626180+00 
0. 244675D+00 
0.2565110+00 
0.2410760+00 
0.2495550+00 
0. 233374D+00 
0 . 2414880+00 
0.2344900+00 
0.2376240+00 
0 . 232665D+00 
0 . 224887D+00 
0.2264590+00 
0.21397BD+00 
0. 1993300+00 
0. 183785D+00 
0. 1677960+00 
0. 15380 1D+00 
0. 132095D+00 


0.437391D-01 
0.756058D-01 
0. 1279440+00 
0. 176326D+00 
0.2140670+00 
0.2309670+00 
0. 247712D+00 
0 . 256568D+00 
0.2608610+00 
0. 263076D+00 
0.2630630+00 
0 . 261599D+00 
0 . 258777D+00 
0. 257403D+00 
0 . 25592 1D+00 
0 . 252464D+00 
0. 250287D+00 
0. 2460B7D+00 
0 . 245 142D+00 
0.2410950+00 
0. 239155D+00 
0. 233013D+00 
0. 227697D+00 
0. 214610D+00 
0.2131860+00 
0. 196894D+00 
0. 180655D+00 
0. 162783D+00 
0. 146695D+00 
0. 1267120+00 


0.2591940-02 
0.4322160-02 
- . 3 15855D-02 
-.8383760-02 
-  .  174805D-01 
-.5459170-02 
- . 90B774D-02 
-.4815630-02 
0. 156848D-02 
0. 178068D-01 
0. 136885D-01 
0.5121130-02 
0 . 384 1210-02 
-.  1272830-01 
0.5897820-03 
-.  1138790-01 
-.7321910-03 
-.  127 13  ID-01 
-.3654260-02 
- .6605450-02 
-.  1530930-02 
- . 347925D-03 
-.2810290-02 
0. 1184870-01 
0.7922930-03 
0.2435650-02 
0 . 3 12989D-02 
0 . 50 1307D-02 
0 . 7 10610D-02 
0 . 53829  ID-02 


•  v< 


m 


x  benzene 


0 . 1604020+00 
0. 1220460+00 
0. 9092800-01 
0.5931000-01 
0. 1949200-01 

The  standard  0 


mo  1  nr  Ve 


Ve  fit 


Error 


0. 1063050+00 
0.890580D-01 
0.6741800-01 
0 . 3863400-01 
0.9062000-02 


0.  1031 10D+00 
0 . 795076D-01 
0 . 59B027D-0 1 
0 . 3934530-01 
0 . 1305400-01 


0.3194770-02 
0.955043D-02 
0 . 761532D-02 
-.71 1347D-03 
- .3992020-02 


viation  for  this  set  of  fitting  constants  is:  0.801066D-02 


The  condition  number  for  the  A  matrix  is: 


13.978 


gasas 


ft 


How  many  curve  fitting  constants  do  you  want? 
(must  De  between  2  and  5) 

2 


■VV  V "J 

,V 

■  vv 


wnat  is  the  name  of  the  data  file  to  be  used? 
bztce25 

(Data  must  be  stored  in  x.y  pairs  with  format  f 7 . 6  .  lx , f 7 . 6 . 
The  first  line  of  the  file  must  contain  the  number 
of  data  points  in  the  file,  with  format  13) 


gamma  1  = 
gamma  2  = 

x  benzene 

0. 130130D-C1 
0 .3754400-0  1 
0.5521000-01 
0 . 78  11400-01 
0 . 9904500-01 
0  .  1 182960*00 
0 .  1472990  +  00 
0 .  1776410+00 
0 .  193  1070  +  00 
0 . 2208740+00 
0 . 2392330+00 
0 . 2588030+00 
0 . 28  17500  +  00 
0 . 3  16878D  +  00 
0.3313990+00 
0.3533900+00 
0.3651730+00 
0 . 3823140+00 
0.3957710+00 
0.4105790+00 
0 .4247610+00 
0 . 4344730+00 
0 .4393690+00 
0 . 4454340+00 
0 . 4492270+00 
0 . 4569270+00 
0.4571390+00 
0 . 4650950+00 
0 .4713760+00 
0 .4767920*00 
0.4881420+00 
0 . 4889200+00 
0 . 4953540+00 
0.4979260+00 
0.5016450+00 
0.5122410+00 
0.5123820+00 
Q .5236970+00 
0 . 5331660  +  00 


0.87624 
-0 . 23889 


molar  Ve 

0 . 0000000+00 
0. 2107400-01 
0 . 3803400-01 
0 . 47  18000-0  1 
0.6102300-01 
0.6805400-01 
0 . 9083700-01 
0 . 1  18  1680  +  00 
0  1  1480  10  +  00 
0.  1315960  +  00 
0. 1308950+00 
0 .  1572550  +  00 
0 . 1588380  +  00 
0 . 1775590+00 
0 . 1782000+00 
0 .  185  2 19D  +  00 
0 . 1922480+00 
0 . 1976480+00 
0 . 2088610+00 
0 . 2095430+00 
0 . 2084020+00 
0.21141 10  +  00 
0 . 2091230+00 
0 . 2228480+00 
0 . 2 117620+00 
0 . 2 188090+00 
0 . 2120860+00 
0 . 2 187550+00 
0 . 227651 D+00 
0 . 2224080+00 
0 . 2241920+00 
0 . 22169  10  +  00 
0 . 2255060+00 
0 . 2299540+00 
0 .2178960+00 
0 . 2221750+00 
0 . 2270930+00 
0 .2251100+00 
0 . 2292540  +  00 


0 . 826384D-02 
0 . 2367840-01 
0 . 3462 12D-0  1 
0.4858450-01 
0.6109660-01 
0 . 7237  19D-0 1 
0.888917D-01 
0. 105506D+00 
0 . 1  13686D  +  00 
0. 1278410  +  00 
0.  1368010+00 
0. 1459780+00 
0. 15622DD+0D 
0. 1 7  J  5  190+00 
0. 176303D+00 
0. 1842190  +  00 
0. 1881980+00 
0 . 193646D  +  00 
0. 1976320+00 
0. 2017 1 4D+00 
0. 2053160+00 
0 . 2076050+00 
0 . 208703D+00 
0 . 2 1001 1D+00 
0 . 2 107990+00 
0 . 2 1 2328D+00 
0 . 2 12368D+00 
0 . 2 13843D+00 
0 . 214934D+  00 
0 . 2 158  2  2D+  00 
0. 217521D+00 
0 . 21763QD+Q0 
0 . 2184860+00 
0 .2188080+00 
0 . 219254D+00 
0 . 2203900+00 
0 . 2204040+00 
0.2213920+00 
0 . 2220400+00 


- .8263840-02 
- . 2604400-02 
0.3412790-02 
- . 140445D-02 
- . 7360560-04 
- .431786D-02 
0 . 194530D-02 
0 . 126624D-01 
0 . 1 1 1536D-02 
0 . 3755090-02 
-.5905680-02 
0.  1127690-01 
0. 26180BD-02 
0 . 6039640-02 
0 .  1897  100-02 
0.9996710-03 
0 . 4050 1  ID-02 
0 . 4002280-02 
0. 1122830-01 
0 . 7828880-02 
0. 308582D-02 
0 . 3805910-02 
0 . 4 198 12D-03 
0 . 1283700-01 
0 . 962966D-03 
0.6481400-02 
- . 2823210-03 
0 . 491 153D-02 
0. 1271680-01 
0 . 6586100-02 
0 . 6670820-02 
0.4061400-02 
0 . 7019830-02 
0 . 1 114550-01 
-  .  1358100-02 
0 . 1785050-02 
0.6689290-02 
0.3717920-02 
0. 7213760-02 
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How  many  curve  fitting  constants  Oo  you  want? 
(must  be  between  2  and  5) 

4 


What  is  the  name  of  the  d.'.ta  file  to  be  used? 


(Data  must  be  stored  in  x,y  pairs  with  format  f  7 . 6 , 1  x  ,  f  7 . 6  . 
The  first  line  of  the  file  must  contain  the  number 
of  data  points  in  the  file,  with  format  13) 


bz  t  ce30 


gamma  1  = 
gamma  2  = 
gamma  3  - 
gamma  4  = 


0.89426 
-0.30734 
0. 10541 
0 . 13159 


x  benzene 


molar  Ve 


Ve  fit 


Error 


0 . 9894600+  00 
0 . 9823640+00 
0. 9724050+00 
0.9666690+00 
0.9406910+00 
0.9380090+00 
0.9229290+00 
0.9003600+00 
0.8776520+00 
0.8761500+00 
0.8602730+00 
0 . 844404D+00 
0.8229840+00 
0 . 80  1236D+00 
0 . 784230D+00 
0 . 7574690+00 
0 . 70934 1D+00 
0.6998650*00 
0.6784220+00 
0.6649030+00 
0 . 639829D+00 
0.6228580+00 
0.6076710+00 
0.5917080+00 
0.5769180+00 
0.5596020+00 
0.5521870+00 
0.5432400+00 
0 . 54 165 1D  +  00 
0 . 533647D+00 
0.5299170+00 
0.5203550+00 
0.5148790+00 
0.5105610+00 
0 . 4958420+00 


0 . 1124300-01 
0. 1810700-01 
0.3955600-01 
0.2984000-01 
0.6863800-01 
0.6566100-01 
0.9166800-01 
0. 1028900+00 
0 . 1254200+00 
0. 1218290+00 
0 . 1372650  +  00 
0. 1461670+00 
0. 1564120+00 
0. 1725570+00 
0. 1828300+00 
0. 192822D+00 
0.2078770+00 
0.2101630+00 
0.2236310+00 
0.2218160+00 
0.2296470+00 
0.2259170+00 
0 .2273780+00 
0 . 228442D+00 
0.2274000+00 
0.2273200+00 
0.2301630+00 
0 . 2263530+00 
0.2302650+00 
0. 220952D+00 
0 . 23 1456D+00 
0.2253740+00 
0 . 2276050+00 
0.2221430+00 
0.2207750+00 


0. 122299D-01 
0.202826D-01 
0.3133500-01 
0.375670D-01 
0.6454710-01 
0.6721500-01 
0.8180140-01 
0. 102306D+00 
0. 121320D+0Q 
0. 122520D+00 
0. 134774D+00 
0. 1462280+00 
0. 1604400+00 
0. 1734160+00 
0. 182558D+00 
0 . 195 193D+00 
0.2127240+00 
0. 215421D+00 
0.2206480+00 
0 . 2  23335D+00 
0. 2271240+00 
0 . 228843D+00 
0.2298290+00 
0. 2303230+00 
0. 230307D+00 
0. 22973 1D+00 
0.2293090+00 
0.2286630+00 
0 . 228549D+00 
0.2278090+00 
0 . 227433D+00 
0 . 226359D+00 
0. 225674D+00 
0 . 225099D+00 
0. 2229120+00 


-.9869340-03 
- . 2 17564D-02 
0.822 104D-02 
- . 772700D-02 
0 . 409092D-02 
-. 135402D-02 
0.986656D-02 
0 . 5B4055D-03 
0 . 409986D-02 
- . 69 1479D-03 
0 . 249 1 23D-02 
- . 612595D-04 
- . 4028 15D-02 

-  .  859403D-03 
0. 272153D-03 
-.2370560-02 
- . 484723D-02 
-.5258270-02 
0.2983470  02 

-  .  151867i  J 2 
0. 252302D-02 

-  .  292620D-02 
-.2450510-02 

-  .  188 125D-02 
-.  290651D-02 
- . 24 1 139D-02 
0.8542810-03 
- . 230954D-02 
0. 17 1589D-02 
- . 685657D-02 
0.402348D-02 
-.9846710-03 
0. 193  107D-02 

-  .  29559  ID-02 
- . 213748D-02 


Ill 


*  Benzene  molar  Ve  Ve  fit 


Error 


0.495406D+00 
0 . 4877940+00 
0.4828 130+00 
0 . 4776850+00 
0.4727170*00 
0.4675540*00 
0.4580930*00 
0.4522000*00 
0.4453140*00 
0.4394220*00 
0.4345690+00 
0 . 4257420+00 
0 .4143170+00 
0 . 3990960+00 
0.388208D+00 
0 . 367306D+00 
0 . 348663D+00 
0 . 326372D+00 
0.3106190+00 
0. 284206D+00 
0.2612640+00 
0 . 2407270*00 
0.2154950+00 
0. 192976D+  00 
0. 164728D+00 
0. 1140820+00 
0.7405700-01 
0.3991800-01 
0. 1822800-01 


0.2238100+00 
0 . 2254530+00 
0 . 2224690+00 
0.2170900+00 
0. 2230460+00 
0.2181410+00 
0.2168510+00 
0 .2167560+00 
0. 2111190+00 
0 .2122240+00 
0.207914D+00 
0 . 209205D+00 
0.2056340+00 
0.2010640+00 
0 . 1979120+00 
0 . 1923940+00 
0. 1857620+00 
0. 179547D+00 
0. 1685350+00 
0. 1634100+00 
0. 1522710*00 
0. 1429610*00 
0. 1289060+00 
0. 1193860+00 
0. 1114010+00 
0.7620200-01 
0.5877700-01 
0 . 2975200-01 
0 . 140990D-01 


0 . 222842D+00 
0. 221573D+00 
0.2206950+00 
0. 219753D+00 
0 . 218803D+00 
0 . 217  778 D+00 
0 .2158030+00 
0. 2145110+00 
0 . 212942D+00 
0.21 1550D+00 
0 . 210370D+00 
0 . 208 148D+00 
0.20513 1D+00 
0.2008740+00 
0 . 197670D+00 
0. 1911610+00 
0. 1849800+00 
0. 177 148D+00 
0. 17  1334D+00 
0. 161084D+00 
0. 151677D+00 
0. 142855D+00 
0 . 1314840+00 
0 . 1208120+00 
0. 1066680+00 
0 . 788673D-0 1 
0 .541921D-01 
0 . 307B34D-0 1 
0 . 1456180-01 


0 . 967507D-03 
0. 387963D-02 
0. 177372D-02 
-.2662690-02 
0.4245190-02 
0 . 362953D-03 
0. 1048060-02 
0.2245420-02 
-.  1822570-02 
0.6742600-03 
-  .  2455880-02 
0. 105690D-02 
0 . 502708D-03 
0. 189559D-03 
0. 242495D-03 
0. 1232940-02 
0 . 78 19B4D-03 
0. 239902D-02 
- . 279929D-02 
0. 232599D-02 
0 . 594324D-03 
0. 1055650-03 
- . 257B25D-02 
-. 1426420-02 
0.47331BD-02 
- . 266533D-02 
0 . 45B487D-02 
-. 1031360-02 
-.4628260-03 


The  standard  deviation  for  this  set  of  fitting  constants  is:  0.324163D-02 


The  condition  number  for  the  A  matrix  is:  102.370 


112 
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Mow  many  curve  fitting  constants  Oo  you  want? 
(must  be  between  2  and  5) 

5 


What  iw  the  name  of  the  data  file  to  be  used? 


(Data  must  be  stored  in  x.y  pairs  with  format  f 7 . 6 , lx , f 7 . 6 . 
The  ’:rst  line  of  the  file  must  contain  the  number 
of  da*a  points  In  the  file,  with  format  13) 


gamma  1  =  0.85454 
gamma  2  =  -0.33061 


gamma  3  =  -0.05882 
gamma  4  =  0.00471 


gamma  =  = 


x  ter 2ene 


molar  Ve 


0 . 96Sl;  ^9D+00 
0.9514720+00 
0. 9184370+00 
0.890; 20D+00 
0 . 869-680+00 
0.8422  20D+00 

0 .8237350+00 
0.7925270+00 
0 . 765307D+00 
0.738845D+00 
o.7oe:-34D+oo 
0 . 679- 24D+00 
0 . 6558  46D+00 
0 . 6336010+00 
0.6131 55D+00 
0.5874950+00 
0.5595440+00 
0.534600D+00 
0 . 5293760+00 
0 . 5 1 34 66D+00 
0.5124690+00 
0 . 493974D+00 
0.4935670+00 
0.4792=00+00 
0.4677090+00 
0.4605520+00 
0.446-610+00 
0.4423280+00 
0.4187310+00 
0.4143600+00 
0.3925270+00 
0.3906180+00 
0.3654690+00 
0.3644790+00 


0.4341600-01 
0.6482000-01 
0.989470D-01 
0.1199820+00 
0.1354360+00 
0.1519760+00 
0. 1651900+00 
0.1827080+00 
0. 1932520+00 
0. 1973 15D+00 
0 . 208375D+00 
0 . 2 1 1747D+00 
0.2174830+00 
0.2125330+00 
0. 2193820+00 
0.2216910+00 
0. 217339D+00 
0 . 2 17635D+00 
0.2147240+00 
0. 214054D+00 
0. 211640D+00 
0. 211692D+00 
0. 212702D+00 
0.2050730+00 
0. 204068D+00 
0. 202780D+00 
0.2038160+00 
0. 1926700+00 
0. 197095D+00 
0. 1919830+00 
0. 1891770+00 
0. 190993D+00 
0. 1796080+00 
0 . 1807060+00 


0 . 432678D-0 1 
0 . 65 1B25D-01 
0.9858640-01 
0. 12 1693D+00 
0. 136899D+D0 
0.  153309D+00 
0.  163292D+D0 
0. 177835D+00 
0. 1885340+00 
0. 197356D+00 
0 . 205736D+00 
0.2119890+00 
0. 215810D+00 
0. 218445D+00 
0. 220002D+00 
0. 220776D+00 
0. 220123D+00 
0. 218235D+00 
0 . 217686D+00 
0. 215691D+00 
0 . 215555D+00 
0 .2126050+00 
0. 212533D+00 
0 . 2098 17D+00 
0. 207368D+00 
0 . 20573BD+00 
0. 202234D+00 
0.2012 140+00 
0. 1946 17D+00 
0. 1933 15D+00 
0. 1864 19D+00 
0.1858560+00 
0 . 1771400+00 
0. 1767870+00 


0.1482420-03 
- . 3624750-03 
0.360564D-03 
-.  1711330-02 
-.  146329D-02 

-  .  133307D-02 
0. 189772D-02 
0.4B7340D-02 
0.4717640-02 
- . 40B940D-04 
0. 263893D-02 
-.2421690-03 
0.  1672980-02 
- . 59 1 246D-02 
- . 620 193D-C3 
0 . 9 15049D-03 
- . 278402D-02 
-.6001810-03 
- . 2962220-02 
- . 1637230-02 
- . 39 1460D-02 
-.9126640-03 
0 . 1692140-03 
-.4743850-02 
- .  330042D-02 
-.2957750-02 
0. 158 180D-02 

-  .  854402D-02 
0.2478460-02 
-. 133248D-02 
0. 275795D-02 
0.51372  ID-02 
0. 246812D-02 
0 . 39  192  ID-02 
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fit 

TO !• 


This  program  is  designed  to  fit  the  data  found  in  a  user  specified 
data  file  to  a  smoothing  equation  which  is  standard  for  excess 
volume  data.  The  data  file  must  be  set  up  in  a  specific  format, 
which  may  be  accomplished  by  using  the  supp 1 i men t a ry  program 
titled  "makefile"  to  construct  the  data  set.  Actually,  the 
program  fits  the  data  for  Ve/(xlx2)  and  then  compensates  after 
fitting  to  arrive  at  the  final  Ve-fit  data  and  the  correct 
standard  deviation.  This  method  was  chosen  to  lessen  the  severity 
of  the  i 1 1 -cond 1 t 1 oni ng  of  the  matrix  while  performing  the  least- 
squares  analysis  on  the  data.  The  program  utilizes  IMSL  subroutines 
to  calculate  the  condition  number  of  the  "A"  matrix  and  to  solve  the 
linear  system  of  equations.  The  user  must  specify  how  many  fitting 
constants  are  desired  in  the  analysis.  (May  be  between  2  and  5). 

Variable  dictionary: 

a:  coefficient  matrix  from  |a|x=P 

b:  solution  vector  from  |A|x=b 

cond:  the  LI  type  condition  number  of  matrix  A 

error:  difference  between  the  actual  and  smoothed  excess  volume 

fac:  matrix  of  the  Lu  f ac tor i za t i on  of  the  A  matrix.  Not  used, 

gamma:  vector  of  I eas t -squares  fitting  constants, 

infoflle:  user  specified  data  file  containing  experimental  data 

ipvt:  vector  containing  pivoting  information.  Not  used  here, 

k:  user  specified  number  of  fitting  constants 

n:  number  of  data  points  in  the  data  file 

rcond:  inverse  of  the  LI  condition  number 

sse :  sum  of  the  square  error  for  all  data  points 

v:  experimental  excess  volume  data 

vefit:  excess  volume  from  the  smoothing  equation 

wkarea:  the  work  area  required  by  IMSL  version  9  routine  LEUT2P 

X:  mole  fraction  of  more  volatile  component 

*  Dimension  variables  and  declare  double  precision 

real *8  a(5,5),fac(5,5) . x(200) , v( 200) ,b(5) .gamma (5) , vef i t ( 200) 
rea I  *  8  error(200) ,wkarea( 100) , sse , rcona , cond 
integer*4  ipvt(5),k,n 
character*14  infofile 

*  Fill  data  vectors  and  arrays  with  zeros 

data  a/ 25*0 . 0d0/ 
data  x / 200*0 . 0d0/ 
data  v/ 200*0 . OdO / 
data  b/5*0.0d0/ 
data  gamma/5*0 . 0d0/ 
data  vef i t/200*0 . 0d0/ 
data  error/200*0 . 0d0/ 

*  Interact  with  user  for  number  of  fitting  constants  and 

*  name  of  data  file 

print  *,  'How  many  curve  fitting  constants  do  you  want?' 
print  *,  '(must  be  between  2  and  5)' 
print  * 


print  *,  'What  is  the  name  of  the  data  file  to  be  used?' 
print  • 

print  *,  '(Data  must  be  stored  in  x.y  pairs  with  format  f7.6.1 
print  *,  '  The  first  line  of  the  file  must  contain  the  number' 
print  *,  '  of  data  points  in  the  file,  with  format  13)' 
print  * 

read  '(A)',  infofile 
print  * 
print  * 


Fill  x  and  V  files  with  experimental  data 

open(unit=  1,  f i 1 e=  infofile,  status=  'old') 
read  (  1  ,  10) , n 
format (  1  x , i 3 ) 

read  (1,20),  (  x  (  i  )  ,  v  (  i  )  ,  i  =  1  ,  n ) 

f  o  rma  t(lx,f7.6,lx,f7.6) 
c I ose ( un i t  =  1 ) 

Convert  Ve  data  to  form  Ve/(xlx2) 
do  50  i = 1 . n 

v(i)=v(i)/(x(1)-x(i)**2) 

cant i nue 

Calculate  coefficient  matrix 

do  51  j  =  1 , k 

do  52  1=1. k 

do  53  i = 1 , n 

a(j,l)=a(j,l)+(l-2*x(i))**(j+l-2) 

conti nue 

cont i nue 

cont i nue 

Calculate  solution  vector 

do  54  1  =  1  , k 

do  55  i  =  1 ,  n 

b( 1 )=b( 1 )*v(i )*( X-2*x( i ) )*»( 1-1) 

gamma( 1 )=b( I ) 

conti nue 

cont i nue 

Call  IMSL  subrout i nes ...  For  convenience,  both  version  9 
and  version  10  function  calls  are  included 

call  DLFCRG  (K , A , 5 , FAC ,5 , IPVT ,RC0ND) 
call  OLSARG  ( k  ,  a  ,  5  ,  b  ,  1  , gamma ) 

C0ND=1 . 0D0/RC0ND 
(IMSL  version  10) 


*  call  LEQT2F ( a , 1 . k , 5 , gamma , 4 , wkarea . i er ) 

*  (IMSL  version  9) 

*  Print  least-squares  fitting  constants 

print  104,  ( i , gamma ( i ). i = 1 , k ) 

104  f ormat ( lx gamma  ',11.'  =  ',f9.5) 

print  • 
print  * 

*  Calculate  Ve  fit  for  each  data  point  and  correct  to  the  form 

*  Ve  from  the  form  Ve/(xlx2).  Also  calculate  error  for  each 

*  data  point  and  accumulate  square  error. 

do  204  i  =  1  ,  n 

do  205  j -  1 . k 

vef i t ( i )=vef i t ( i )+ gamma (j)*(l-2*x(i))»*(j-l) 

205  continue 

vefit(i)=vefit{i),(x(i)-x(i)**2) 
v ( i  )=v( i  )*(x(  i  )-x( i )*«2) 
error(i)=v(i)-vefit(i) 
sse=sse»error(  i  )  *  *  2 

204  conti nue 

*  Calculate  standard  deviation  for  the  data  set 

s igma= (sse/(n-k) .5 

»  Print  data  summary  point  by  point 

print  *,  '  x  benzene  molar  Ve  Ve  tit 

cError ' 

print  *,' _ ' 

print* 

print  106,  ( x ( i  ),  v ( i ), vef i t ( i ), error ( i ),  i  =  1 , n) 

106  f ormat ( lx . 4 (d 12 . 6 , 4x ) ) 

print  • 
print  * 
print  * 

*  Print  standard  deviation  information 

print  108 ,  si gma 

108  f ormat ( lx ,' The  standard  deviation  for  this  set  of  fitting  constant 

cs  is:  ' , d 12 . 6 ) 

print  ♦ 
print  * 

*  Print  the  condition  numoer  of  the  A  matrix 

PRINT  109.  C0N0 

109  FORMAT ( IX ,' The  condition  number  for  the  A  matrix  is:  ',f8.3) 
end 


This  program  is  designed  to  supplement  the  program  named  curvf  it. 

Its  purpose  is  to  create  a  data  file  of  experimental  data  points 
that  is  in  a  proper  format  to  be  accessed  by  the  curvf i t  program. 

This  program  interacts  with  the  user  to  determine  the  name  of  the 
data  file,  and  then  prompts  the  user  to  enter  the  data  in  the 
proper  format.  The  data  are  arranged  here  as  a  text  file,  but  are 
read  by  the  program  curvf it  as  numeric. 

Variable  dictionary: 

infofile  :  the  user  specified  name  of  the  data  file  being  made 
k  :  the  number  of  data  points  to  be  included  in  the  file 
line  :  the  information  from  a  single  data  point. 

Declare  variables  and  dimension  them 

character  *20  line(200).  infofile 

Ask  the  user  for  the  name  of  the  data  file  to  be  created 

print  *,  'What  is  the  name  of  the  file  you  want  to  create?' 
read  '(A)',  infofile 

Read  in  the  data  elements  to  the  file,  number  of  data  points  first. 

open  (unit  =  1,  file  =  infofile,  status  =  'new') 

print  ».  'How  many  data  points  will  be  in  this  file?' 

read  *  ,  k 

write  (1,10),  k 

format (  1  x , i 3 ) 

Prompt  the  user  to  enter  the  data  points  in  the  correct  format. 

print  *,  'Enter  the  data  points  as  x,y  pairs,  where  x  is' 
print  *,  'the  mole  fraction  of  the  more  volatile  component' 
print  *,  'and  y  is  the  experimental  molar  excess  volume.' 
print  * 

print  *,  '(There  should  be  a  space  before  the  x  points  and' 
print  *,  '  a  space  between  the  x  and  y  points,  and  their  formats 

print  *,  '  should  be  f7.6.)' 

print  ♦ 

print*,'  x  Ve ' 

pr i nt  * 

Read  in  the  data  points. 

do  100  i  =  1 ,  k 

read  20 ,  I i ne  (  i  ) 
write  (1,20)  line(i) 
format ( lx , a20 ) 

0  cont i nue 


c I ose  (unit  =  1 ) 


APPENDIX  E 


BURET  CALIBRATION  DATA 


TABLE  31 


EXPERIMENTAL  DATA  FOR  BURET  CALIBRATION 


RUN  1 


POINT  NUMBER  WT  OF  HG 

ADDED 


HT  OF  HG 
IN  BURET 
(  CM  > 


REF  LINE 
HEIGHT 
<  CM  ) 


INITIAL 

1 

2 

3 

4 

5 


136. 92488 
25. 31131 
20. 76882 
17. 22815 
22. 66185 
15. 37205 


0. 2723 
2. 0117 
3. 9717 
5. 5975 
7. 7166 
9.  1562 


O. 1575 
-O. 5710 
-O. 5063 
-0. 5184 
-O. 5101 
-O. 5276 


RUN  2 


POINT  NUMBER 

WT  OF  HG 

HT  OF  HG 

REF  LINE 

ADDED 

IN  BURET 

HEIGHT 

(  G  > 

(  CM  ) 

(  CM  > 

INITIAL 

147. 70806 

0. 6726 

-0. 4971 

1 

14. 28067 

1 . 9920 

-0. 5126 

2 

14. 22029 

3. 3230 

-0. 5250 

3 

12. 84957 

4. 5506 

-0. 5036 

4 

13. 54833 

5.  7898 

-0. 5378 

5 

13. 78516 

7. 1099 

-0. 5104 

6 

13. 56552 

8. 3903 

-0. 5119 

7 

17. 10570 

10. 0010 

-0. 5194 

RUN  3 


POINT  NUMBER 

WT  OF  HG 

HT  OF  HG 

REF  LINE 

ADDED 

IN  BURET 

HEIGHT 

(  G  ) 

<  CM  ) 

(  CM  > 

INITIAL 

162. 67139 

2. 0926 

-0. 4613 

1 

12. 03903 

3. 1739 

-0. 5217 

2 

13. 28431 

4. 4089 

-0. 5414 

3 

12. 06771 

5. 5758 

-0. 5109 

4 

14. 32788 

6. 8920 

-0. 5266 

5 

15. 23696 

8. 3225 

-0. 5400 

6 

15. 02251 

9. 7895 

-0. 4939 

APPENDIX  E-2 


CALCULATED  DATA 
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Cl 

ca 

[  <  Cl  . 

i 

-  CR±  ) 

-  (C^  -  CRj^ 

)  3 

C2 

ca 

[  (C2± 

-  CR±  ) 

-  <C2jl  -  CR1 

)  ] 

C 

ca 

t<Ci  - 

CR± )  - 

<C1  -  CR^ ) 3 

Where:  V  represents  a  volume 

B  refers  to  buret  measurements 

BR  refers  to  buret  reference  line  measurements 

Cl  refers  to  measurements  on  capillary  Cl 

C2  refers  to  measurements  on  capillary  C2 

C  refers  to  measurements  on  capillary  C 

CR  refers  to  capillary  reference  line 
measurements 


subscript  i  indicates  an  initial  measurement 
subscript  1  indicates  a  load  measurement 
subscript  ca  indicates  cross  sectional  area 
subscript  BB  indicates  buret  bulb 
subscript  sc  indicates  solvent  charged 

For  run  #17 : 

V  =  10.0069  «■  .78452  (7.53-7.41) 

£C  ♦  .01773  [(.  6344-1.  4803) - (4. 8066-1. 6693) 3 

♦  .01776  [(. 5530-1. 4803) -< 4. 7747-1. 6693) 3 
-  .01767  C <3. 6977-1. 4803) - <4. 8263-1. 6692) 3 
=  9. 9422  cc 

2.  holes  of  solvent  charged  in  the  mixing  bowl: 

n  =  V  /  v 
sc  sc  solvent 

where:  n  is  the  number  of  moles  of  solvent  charged 

SC 

v  .is  the  specific  volume  of  the  solvent 

solvent 

n  =  9.9422  /  87.8134  =  .1132  gmole 
sc 

The  rest  of  the  data  will  be  illustrated  by  Figures  17,  18, 

and  19. 

3.  Change  in  height  of  mercury  in  buret,  B: 

del  B  =  <B2  -  BR2)  -  (B1  -  BR^ ) 

del  B  =  (8.16  -  7.42)  -  <7.81  -  7.42)  =  .35  cm 

4.  Change  in  height  of  mercury  in  capillary  Cl: 

del  Cl  =  <C12  -  CR2)  -  ( Cl :  -CR1) 

del  Cl  =(3.6384  -  1 . 4938 ) - < 3.  671 1  -  1.4998)  =  -.0267  cm 


FIGURE  18 


RUN  #17,  DATA  POINT  #2  BEFORE  PRESSURE  CORRECTION 


*  nirn¥ru 


FIGURE  19 


RUN  #17,  DATA  POINT  #2  AFTER  PRESSURE  CORRECTION 


C  =  4.1257 


Ux  =  27.65 


U  =  51.60 


m 


pi 

Kw 

m 


-v.;- ,, 

W 

w. 


ii 


p;: 

mb 


M& 


Change  in  height  of  mercury  in  capillary  C 2: 
del  C 2  =  <C22  -  CR2>  -  <  C21  -  CR  1 ) 

del  C2  =  (3.6239  -  1.4936)  -  <3.6494  -  1.4996) 

=  - . 0195  cm 

Change  in  height  of  mercury  in  capillary  C: 

del  C  =  <C„„  -  CR_ )  -  <C„,  -  CR, ) 

B2  2  B1  1 

del  C  =  <4.2256  -  1.4936)  -  <3.9899  -  1.4996) 

=  .2417  cm 

Change  of  height  of  mercury  in  mixing  bowl: 
del  A  =  (A^  -  CR2)  -  (A  -  CR1  ) 

del  A  =  <3.1104  -  1.4938)  -  <3.1422  -  1.4998) 

=  - . 0258  cm 

Volume  of  solute  added  to  solvent: 

del  V  =  (del  B)*B  *  (del  C1)*C1  *  (del  C2)*C2 

c  ca  ca  c 

del  V  =  (.35)*. 78452  *  ( - . 0267 ) * . 01 773 
c 

*  <-. 0195) *. 01776  =  .27376  cc 

Number  of  moles  of  solute  added  to  the  solvent: 

del  n  =  del  V  /  v 

c  c  solute 

del  n  =  .27376  /  88.627319  =  .003089  gmole 

c 

Total  number  of  moles  of  solute  in  the  mixing  bowl : 

n  TOT_  =  n  TOT,  +  del  n 
c  2  c  1  c 

n  TOT_  =  .003610  +  .003089  =  .006699  gmole 

C  2 

Mole  fraction  of  benzene  in  mixing  bowl: 

x  (Benz)  =  n  /  (n  n  TOT) 

sc  sc  c 

x ( Benz )  =  .1132  /  (.1132  *  .006699)  =  .9441 
Cumulative  height  change  of  mercury  in  capillary  C: 
CT0T2  =  CTOT^  *  del  C 


CTOT_ 


27268  +  .2417 


5144  cm 


1  ^ 


Pressure  corrected  change  in  height  of  mercury  in 


capillary  C: 


del  C  pi  '  CA2  "  CB2 

del  C  pi  =  4.1257  -  4.2256  =  -.0999  cm 


14.  Back-pressure  applied: 

del  U  =  (U,w  -  U  w)  «■  <U  -  U,  > 
lb  rb  ra  la 

del  U  =  (39.70  -  39.75)  +  <51.60  -  27.65)  =  23.90  cm 

15.  Excess  volume  (pressure  corrected):  (Kumaran,  p.  265) 

VE  =  CTOT  *  C  *  Cl  -  (del  C  pi  /  del  U)3 
ca 

*  Cl  -  (del  A  /  CTOT)] 


V  =  .5144  *  .01767  *  Cl  -  (-.0999/23.90)3 

«  Cl  -  ( - . 0258/. 5144 ) 3 
VE  =  . 009585  cc 

16.  Molar  excess  volume  (pressure  corrected): 

V  -  VE  /  <n  *■  n  TOT) 

em  sc  c 

V  =  .009585  /  (.1132  ♦  .006699)  =  .0799  cc/gmole 


B.  Buret  Calibration  Data 

Calibration  was  conducted  at  24. 947  °C.  The  mercury 
density  at  this  temperature  is  13. 5337  g /cc.  This  value  was 
obtained  from  Figure  20. 

1.  Cross  sectional  area  (using  run  #2  point  1): 

A  =  w/  <^3  *  A  h  ) 

where:  A  is  the  cross  sectional  area 

w  is  the  weight  of  mercury  added 
P  is  the  density  of  mercury 
^h  is  the  change  in  height  of  mercury 

=  (2.5046  -  1.1697)  =  1.3349  cm 

A  =  14.28067  /  (13.5337  •  1.3349)  =  .7904662  cm2 


Ww** 

R**.v 

mS, 

m\ 


Bow 

mm 

Grow** 

EiWjtiv 

RSm* 


Average  cross  sectional  area  and  standard  deviations 


( using  run  #2 ) : 


A  = 

/  n 

vher  e 

•  A 

S  e 

IKo 

average  erase  sectional  area 

A  . 

i 

n 

is 

the 

cross  sectional  area  of  data 

point  i 

is 

the 

number  of  data  points  in  the 

run 

Using 

ths» 

statisti 

cal  functions  on  a  Texas  Inst 

ruments 

TI-55  III  calculator, 

T  =  5.4968473  /  7  =  .7852639  cm2 

(T  =  .0033168 
n 

rr  =  . 0035825 
n-  1 

Note:  The  final  A  value  used  w as  that  calculated  from  the 

13  data  points  in  runs  2  and  3  combined.  Run  1  vas 
discarded  due  to  poor  agreement  with  the  other  runs. 

3.  Volume  of  the  bulb  BB  (from  run  #2): 

V  =  Volume  of  mercury  initially  charged  -  volume  of 

DI3 

mercury  above  BR 


Weight 

Height 

of 

of 

mercury 

mercury 

charged  = 
charged  = 
A  = 

147.70806  g 

1.1697  cm  above  BR 
. 7853  cm0 

Weight 

of 

mercury 

above  BR: 

w  =  A  +  *^h 

w  =  .7853  *  13.5337  *  1.1697  =  12.431589  g 

Weight  of  mercury  below  BR  =  weight  of  mercury  charged 

-  weight  of  mercury  above  BR 

Weight  of  mercury  below  BR  =  147.70806  -  12.431589 

=  135. 27647  g 

V  =  weight  of  mercury  below  BR  /  density  if  mercury 

O  D 

V__  =  135.27647  /  13.5337  =  9.9955  cc 
b  b 

Note:  The  volume  used  for  the  bulb  was  the  average  for  runs 


2  and  3. 


APPENDIX  G 


NOMENCLATURE 


C,  Cl,  C2 
C,  Cl,  C2 

ci 

CR 

CTOT 
del  A 
del  B 

del  C,  Cl , C2 
del  C  pi 

del  n 

c 

del  V 

c 

del  U 
„E 


, 

. 

-'mm 

NOMENCLATURE 

Mixing  bowl  on  the  dllatometer 
Average  cross  sectional  area 
Cross  sectional  area  of  data  point  1 
Least-squares  fitting  constants 
Buret  on  the  dllatometer 
Bulb  at  the  bottom  of  buret 
Buret  cross  sectional  area 
Reference  mark  on  the  buret 
Condition  number 

Capillary  tubes  on  the  dllatometer 
Capillary  cross  sectional  area 
Reference  mark  on  capillary  C 
Cumulative  change  in  height  of  mercury  in  C 

Change  in  height  of  mercury  in  mixing  bowl 

Change  in  height  of  mercury  in  buret 

Change  in  height  of  mercury  in  a  capillary 

Pressure  corrected  change  in  height  of  mercury 
in  C 

Number  of  moles  of  solute  added  to  solvent 
Volume  of  solute  added  to  the  solvent 
Back-pressure  applied  to  C 
Excess  Gibbs  free  energy  of  mixing 
Number  of  data  points  used  to  fit  a  curve 
Number  of  coefficients  in  least-squares  fit 


mm 

m\ 


sp 
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i??>  & 


n 


Number  of  data  points  in  buret  calibration  run 


n  TOT 
c 

n 

sc 

P 

SI,  S2 
T 

Tl,  T2 

Ui 

V. 

v1 

V 

( 

T 
v. 
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eT 
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sc 


solute 


solvent 


x  ( Benz  > 


"i 

A* 

Av 

cr 

n 

0"«. 


n  - 1 


r 


Total  number  of  moles  of  solute  in  mixing  bowl 
Number  of  moles  of  solvent  charged 
Total  pressure 

Vent  tubes  on  the  dilatometer 
Temperature 

Taps  on  the  dilatometer 

Back-pressure  reading  i  on  U  tube 

Volume  of  bulb  BB 

Excess  molar  volume  of  mixing 

Maximum  excess  volume  at  temperature  T 

Partial  molar  volume  of  component  i 

Molar  volume  of  component  i  in  solution 

Molar  volume  of  an  ideal  solution 

Excess  volume  of  mixing  for  data  point  i 

Volume  of  solvent  charged 

Specific  volume  of  the  solute 

Specific  volume  of  the  solvent 

weight  of  mercury  added 

Mole  fraction  of  benzene  in  mixing  bowl 

Mole  fraction  of  component  i  in  solution 

Change  in  height  of  mercury  in  calibration  run 

Molar  volume  change  on  mixing 

Population  standard  deviation 

Standard  deviation  of  a  point 

Density 
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